Topological transitions in two-electron quantum dots in magnetic field 
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We use entanglement to study topological transitions in two-electron axially-symmetric parabolic 
quantum dots under a perpendicular magnetic field. At a specific magnetic field value the dot attains 
a spherical symmetry. The transition from the axial to the spherical symmetry manifests itself as 
a drastic change of the entanglement of the lowest state with zero angular momentum projection. 
While the electrons in such a state are always localized in the plane (x-y) before the transition point, 
after this point they become localized in the vertical direction. 

PACS numbers: 03.67.Bg, 73. 21. La, 73.22.-f 



Nowadays there is a growing interest in using quan- 
tum entanglement measures for study ofquantum corre- 
lations in topologically ordered systems ^ . The analysis 
of these systems is a highly non-trivial task due to the 
absence of the order parameter. The main stream of such 
an analysis is focused on many-particle one-dimensional 
systems (of [2l]). Recently, topological phases were stud- 
ied in finite systems such as two-dimensional electrons at 
very high magnetic field (Laughlin states) [sj and weakly 
interacting two-dimensional rotating Bose-Einstein con- 
densate [J]. It is noteworthy that various topological 
phases are exhibited in quantum dots (QDs) at differ- 
ent strengths of the applied perpendicular magnetic field 
. Two-electron QDs being realictic tractable nontrivial 
systems are, in particular, attractive because their eigen- 
states can be obtained very accurately, or in some cases, 
exactly (cf [1, 01)- The goal objective of the present pa- 
per is to demonstrate that quantum entanglement can be 
used to trace a topological transition of two interacting 
electrons confined in a three-dimensional (3D) quantum 
dot under a magnetic field. 

Our analysis is carried out by means of the exact di- 
agonalization of the Hamiltonian 
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Here k = e^/A-KEQEr and Hgpin — g*fiB{si + S2)-B de- 
scribes the Zeeman term, where fiB — eh/2meC is the 
Bohr magneton. As an example, we will use the ef- 
fective mass m* — 0.067me, the relative dielectric con- 
stant Er — 12 and the effective Lande factor g* = —0.44 
(bulk GaAs values). For the perpendicular magnetic 
field we choose the vector potential with gauge A = 
iB X r = ^B{—y,x,0). The confining potential is ap- 
proximated by a 3D axially-symmetric harmonic oscilla- 
tor /7(r) = m*[ujQ {x^+y'^)+u}'^z^]/2, where huiz and hujo 
are the energy scales of confinement in the z-direction 
and in the xy-plane, respectively. 

By introducing the center of mass (CM) and relative 



coordinates: R — i(ri -|- r2) and ri2 = ri — r2, the 
Hamiltonian ([ij, in agreement with the Kohn theorem 
Q, separates into the CM and relative- motion terms 
H = HcM + Hrci (see details in [7j]). The CM term is 
described by the oscillator Hamiltonian with the mass 
Ai = 2m* and frequencies of the one-particle confin- 
ing potential U. The Hamiltonian for relative motion 
in cylindrical coordinates takes the form 
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where /i = m* /2 is the reduced mass, Iz (— —ihd/dipi2) 
is the projection of angular momentum for relative mo- 
tion and pi2 = (a;i2 +yh)^^'^, V12 = arctan(2/i2/a;i2), 
fi2 = {Pi2 + 212)^^^- Here, ujl = eB/2m*c is the Larmor 
frequency, and the effective lateral confinement frequency 
— {lo\ + '^o)^^^ depends through on the magnetic 
field. 

The total two-electron wave function \['(ri,r2) — 
-(/"(ri, r2)x(o'i, 0-2) is a product of the orbital ^'(ri, r2) and 
spin x(^ij ^2) wave functions. Due to the Kohn theorem, 
the orbital wave function is factorized as a product of the 
CM and the relative motion wave functions 



■'/'(ri,r2) = -001^(11.) V'rcl(l"l2) 



(3) 



According to the Pauli principle, the orbital wave func- 
tion must be symmetric (or, equivalently, V'rei(ri2) must 
be even) for the antisymmetric (singlet: S = Ms = 0) 
spin state, and it must be antisymmetric (V'rci(ri2) must 
be odd) for the symmetric (triplet: S = 1, Ms = 0, ±1) 
spin states. Thus, for the relative motion the parity of 
'0roi(ri2) is a good quantum number as well as the mag- 
netic quantum number m, since Iz is the integral of mo- 
tion. 

The CM eigenfunction is a product of the Fock-Darwin 
state (the eigenstate of electron in an isotropic 2D har- 
monic oscillator potential under a perpendicular mag- 
netic field) in the {X, y)-plane and the oscillator func- 
tion in the Z-direction (both sets for a particle of mass 



2 



A4). In this paper we consider the lowest CM eigenstate 
which has the form ■(/'cm(R-) = '>Pcm i-^y^) ''I^cm(^)^ 
where ^^^(X,^) = y^^^ e'^'(''"+^"^ and VSi(Z) = 

{2u}z/TrY^^e~'^'^ (i.e. zero principle quantum numbers), 
with = m*D,/h and (D^ = m*uJz/h. 

Since the Coulomb interaction mixes the eigenstates 
of non-interacting electrons, the eigenfunctions of the 
Hamiltonian for relative motion ([2]) are expanded in the 
basis of the Fock-Darwin states $n,m (pi2, ^12) and oscil- 
lator functions in the zi2-direction (jjn^ (^12) (for a particle 
of mass i-j,), i.e. 

V'rol(ri2) = cS.^n,m(pi2,9?12)(/'n,(2;i2)- (4) 



The coefhcients Cn'^n^ can be determined by diagonalizing 
the Hamiltonian ([2]) in the same basis. Evidently, in 
numerical analysis the basis is restricted to a finite set 
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be, however, large enough to provide a good convergence 
for the numerical results. Since the function '!/'rci(ri2) has 
a definite parity and the parity of the functions $ri.m (f'n^ 
is (—1)'"+"', the index in the expansion ^ takes 
either even or odd values. 

For non-interacting electrons (k = 0) the eigenfunc- 
tions -^roi are simply the basis functions , and, 
therefore, the ground state is described by the wave func- 
tion tpj-ci — $0,0 </'o- When two interacting electron move 
in the external field created by the confining potential 
and the applied, varying steadily, magnetic field, the 
quantum number m of the ground state (in the form ^) 
evolves from zero to higher values as the magnetic field 
strength increases. It results in the well known singlet- 
triplet transitions [lo| . Namely, for a given m the domi- 
nant term in the expansion Q will be $o,m (f'o all 
are even) and the parity of the ground state is (—1)™, 
which determines the total spin to be = ^[1 — (—1)™]. 
Note that the quantum number Ms associated with the 
spin wave function evolves as following: for even m the 
total spin S ~ and, thus, A/5 = 0; for odd m the 
total spin S = 1 and Ms can be —1, or 1. The Zee- 
man splitting (with g* < 0) will lower the energy of the 
Ms — 1 component of the triplet states, while leaving the 
singlet states unchanged. As a consequence, the ground 
state will be characterized by Ms = S. With the increas- 
ing magnetic field the intervals of the triplet states will 
increase at the cost of the singlet ones, and eventually, 
the singlet ground states will be totally suppressed. The 
increase of the magnetic field leads to the formation of 
a ring and a torus of maximal density in 2D- and 3D- 
densities, respectively (see Fig. 4 in Refflll). 



At the value uj^^^ 



q Y^^ the magnetic field 



gives rise to the spherical symmetry {ujz/^ = 1) (with 
ujz > Wo) in the axially- symmetric two-electron QD 
[l3| . This phenomenon was also recognized in the 
results for many interacting electrons in self- assembled 
QDs In the later case it was interpreted as an ap- 
proximate symmetry that had survived from the non- 



interacting case due to the dominance of the confinement 
energy over a relatively small Coulomb interaction en- 
ergy. However, the symmetry is not approximate but 
exact even for strongly interacting electrons, because the 
radial electron-electron repulsion does not break the rota- 
tional symmetry. A natural question arises how to detect 
such a transition looking on the density distribution only. 
The related question is, if such a transition occurs, what 
are the concomitant structural changes? 

To this end we employ the entanglement measure based 
on the linear entropy of reduced density matrices (cf [IBl) 
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where p "'"' and p^r ^"''' are the single-particle reduced 
density matrices in the orbital and spin spaces, respec- 
tively. This measure is quite popular for the analysis of 
the entanglement of two-fermion systems, in particular, 
two electrons confined in the parabolic potential in the 
absence of the magnetic field . 

The trace Tr[pf ] of the two-electron spin states 
with a definite symmetry xs,Ms has two values: (i) 1/2 
if Ms = (anti-parallel spins of two electrons); (ii) 1 if 
Ms = ±1 (parallel spins). The condition Ms — S — 
i[l „ (-1)™] yields 



Tr[p(^f^")'] = i(l + |A'/5l) = 



3-(-l)^ 



(6) 



The trace of the orbital part Tr[pr°'^'''' ] 
Tr[p(°'^'')'] = ydridr/dr2dr2'V(ri,r2)^*(r/,r2) 

r{v,,vi)^(v;,vi). (7) 

is more involved. Indeed, in virtue of Eqs.Q, (|4]), one 
obtains 
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where 



/(ni, 712, na, 714; m) = j dridr[dr2dr2 
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(here are vectors in the a;?/-plane) and 
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FIG. 1: (a) Entanglement of the lowest state with m = at 
Rw = 2 and various ratios cJz/tJo as functions of the parame- 
ter ujl/ojo- The circles denote the values of ujl/ujo when QDs 
with given ratios oJz/ojq become spherically symmetric, (b) 
The relative strength of the Coulomb interaction R'^^^ /Rw 
(solid line) and Rq^-' /Rw (dash-dotted line) for the lowest 
state with m = at various ratios uJz/oJo as functions of the 
parameter cjl/cjo- 



The magnetic field dependence of the entanglement £ 
naturally occurs via inherent variability of the expansion 
coefficients. The values of the / and J integrals for any 
choice of indices can be determined analytically, vifhich 
simplifies the numerical calculations. 

For our analysis it is convenient to use the so-called 
Wigner parameter Rw — {k/loj/huoo = lo/a* , a measure 
of the Coulomb interaction strength rela tive to th e con- 
finement strength (cf [7]). Here, = ^JhJmFujQ is the 
oscillator length and a* ~ /km* is the effective Bohr 
radius. For our choice of the parameters (GaAs) and 
for the confinement frequency Hujq w 2.8 meV we have 
Rw ~ 2. The numerical analysis demonstrates a good 
convergency for the basis with Umax = n™°-^ — 4. 

In the absence of the magnetic field {B — 0) the en- 
tanglement decreases if the ratio uJz/ujq decreases from 
oo (2D model) to 1 (spherically symmetric 3D model); 
see Fig. [TJa) at wl/wq = 0. This effect could be ex- 
plained by introducing the effective charge kes [H, 
which determines the effective electron-electron interac- 
tion Vc = fceff/pi2 in the QD. In the 3D dot the electrons 
can avoid each other more effectively than in the 2D one. 
Therefore, the Coulomb interaction has a smaller effect 
on the 3D spectrum (the ratio k^s/k ^ 0.5) in contrast 
to the 2D case when kcs/k = 1. Thus, a decreasing of the 
ratio uJz/ujQ yields an analogous effect as the reduction of 
the electron-electron interaction. 

Fig. [ija) shows the entanglement measure £ of the 
lowest angular momentum state to = as a function 
of the magnetic field (the parameter ujl/uq) at a fixed 
value of Rw and for different ratios Wz/wq. In the 2D 
case the entanglement decreases monotonically with the 
increase of the magnetic field. The constant electron- 
electron interaction becomes relatively weaker, since the 
effective lateral confinement (hfl) increases with the mag- 
netic field. If we introduce the characte ristic length of 
the effective confinement In = y^h/m*fl, the parameter 



Rn = la/ci* (which is equal Rw at B = 0) determines 
the relative strength of Coulomb interaction at a given 
effective confinement. Evidently, i?a decreases with the 
increase of the magnetic field B (see Fig. [IJb) , the line 
labelled by '2D'). In the 3D case, however, the entangle- 
ment decreases until lul — <^^l^, when the spherical sym- 
metry occurs. After this point the entanglement starts 
to increase (see Fig. [Tfa)). 

This behaviour can be explained by the influence of 
magnetic field on the effective strength Rq, which is 
twofold here. Indeed, in the 3D case the magnetic field 
affects the effective charge as well as the effective con- 
finement. For the quasi-2D system of electrons (fi <C Wz) 
the effective charge is k'^^^ = (^12^0) (see Eq. (18) 
in Ref. [l^), where Vc = k/^p\^ -\^ is the full 3D 
Coulomb interaction. The mean value {pv^yc) is calcu- 
lated by means of the eigenstates of Uxcx in the approxi- 
mation of non-interacting electrons. Here, the eigenstate 
is i'fl.m 4>n (for explicit expressions see Eqs. (19), (20) in 
Ref. [i3|)- Thus, for the quasi-2D case the parameter 

^n^^ — {m* /h^VlY/"^ ^off^' used as a measure for 

the relative strength of the Coulomb interaction. 

For ujz (very strong magnetic field) the electrons 
are pushed laterally towards the dot's center. The mag- 
netic field, however, does not affect the vertical confine- 
ment. As a consequence the electrons practically can 
move only in the z-direction and the QD becomes a 
quasi-lD system. In this case a measure for the rela- 
tive strength of Coulomb interaction can be defined as 
R^^^^ — {m* /h'^ujz)^^^ k'^^\ where the effective charge 

for a quasi-lD system is k^^^ = (\zi2\Vc)- It can be 
shown that for the lowest state with m = one obtains 

fc(i°Vfc = (i + y^)"'- 

The quantities Rq^^ and Rq^^ for the lowest state 
with m = 0, as functions of the parameter ljl/ujq (in 
the domains < cj^ < oj^j^^ and ojl > respec- 
tively), are shown in Fig. [Ijb) for different ratios uJz/ujq. 

One observes that the effective strength R^^^ decreases 
with the increase of the magnetic field for different ratios 
Wz/wo, similar to the 2D case. The oppositely ordered 
confinement 17^^'^^ (which is not defined for the 2D case) 
increases with ujl and, therefore, the effective strength 
Rq^^ increases as well. In order to match R^^^ = Rq^'^ 

at ujl = ijJl'^ (i-e. when $7 = w^) the strength R^^"^ 
is scaled by the factor 7r/2. Although at this point the 
3D system is far from the 2D model and from the ID 
model and, as a consequence, R^^^ and i?^'^'' do not 
match smoothly, these two functions taken together give 
a qualitative picture how the effective electron-electron 
interaction in a 3D QD changes with the magnetic field. 

To get deep insight into this transition we calculate 
the probability density |'0(ri2)P and potential surfaces 
for various values of the magnetic field (see Fig. [2]). Since 
the symmetry is exact for any strength of the electron- 
electron interaction at the transition point, in order to 
illuminate the effect, we use Rw = 10. For the magnetic 
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FIG. 2: The probability density |'0(ri2)j^ of the lowest m — 
state (top) and the contour plots of the potential surface 
(bottom) for the QD with ujz/ijJo ~ 2 and Rw = 10 shown 
in (pi2, Zi2)-plane for the cases: (a,d) lul/i^o = (ajz/fi — 
2), (b,e) ojl/ujo = ojf'^/oJo = 1.73205 (cj^/!^ = 1) and (c,f) 
ujl/ujo = 2.29129 (w^/fi = 0.8). 



(fi < Uz) the stationary point pi2 = poj -^12 = is the 
minimum of the potential surface (see Fig. djd)). Here 
po = (fc//if2^)^/^ [13]. By increasing the magnetic field 
over the value Bsph > ^z) the stationary point trans- 
forms to the saddle point and two new minima appear, 
divided by a potential barrier (see Fig. [H^f)). In other 
words, for m = a bifurcation of the stationary point lo- 
cated at (po, 0) occurs at the value of magnetic field when 
lol ~ w^^'^ (see Fig. [2jb,e)). In the domain VL > uj^, for 
TO = 0, the minima are located at Z12 = ±zo in the Z12- 
axis (pi2 = 0), where zq = {k/ ^loIY^^ . 

Summarizing, we have shown that the 3D approach 
provides a consistent description of the topological tran- 
sition in two-electron QDs under the magnetic field. The 
entanglement of the lowest state with m = 0, being first 
a decreasing function of the magnetic field, starts to in- 
crease after the transition point with the increase of the 
magnetic field. This behaviour is understood as the tran- 
sition from the lateral to the vertical localization of the 
two-electron probability density for this state in the QD. 



field wl < ^^"^ density maximum is located in the 
(a;i2,yi2)-plane (zi2 = 0, see Fig. [SJa)). For > i^f^, 
however, there are two maxima located symmetrically 
along the zi2-axis (pi2 = 0, see Fig. ^c)). The analysis 
of the behaviour of the stationary point of the poten- 
tial V = ^p(r2^p22 + '^1^12) + k/ri2 as a function of the 
magnetic field provides the explanation. For ujl < w^^^ 
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